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The 2-dimensional U(l) gauge-Higgs model with a topological term is a simple example of a lattice 
field theory where the complex action problem comes from the topological term. We show that the 
model can be exactly rewritten in terms of dual variables, such that the dual partition sum has only 
real and positive contributions. Using suitable algorithms the dual formulation allows for Monte 
Carlo simulations at arbitrary values of the vacuum angle. We demonstrate the feasibility of the 
dual simulation and study the continuum limit, as well as the phase diagram of the system. 


While working on the reply to the referee report the sad news reached us that our friend and 
co-author Michael Muller-Preussker has passed away on 12th of October 2015. Michael has devoted 
50 years of his life to physics as a very successful researcher, but also was a colleague for whom 
service to the community was an important duty. Michael was an inspiration for many and his 
untimely early death is a big loss for all of us. 


I. INTRODUCTORY REMARKS 


Systems with a vacuum term have a complex ac¬ 
tion problem that is similar to the one for systems 
with finite chemical potential. For the latter case it 
was recently shown that some lattice models can be 
rewritten exactly to new degrees of freedom, so-called 
dual variables. In terms of the dual variables, which 
are loops for matter fields and surfaces for gauge fields, 
the partition function has only real and positive con¬ 
tributions and Monte Carlo simulations are possible. 
For examples relevant to this project see M- 

For clarity we stress that in the strict sense the vari¬ 
ant of duality applied here is only the first half of the 
conventional duality transformation, as will also be¬ 
come clear below. On conventional duality transfor¬ 
mations there exists extensive literature - see, e.g., 
the review [5j. We remark, that the second step of 
the conventional duality transformation, i.e., satisfy¬ 
ing the constraints by introducing variables on the 
dual lattice is often not possible for the models con¬ 
sidered in m, which, however, is irrelevant for a 
numerical simulation. 

In these notes we apply the same dualization tech¬ 
niques as in M to a simple system with a topological 
term, the U(l) gauge-Higgs model in two dimensions. 
We show that a dual representation with only real and 
positive contributions is possible for arbitrary values 
of the vacuum angle 8 and implement a suitable Monte 
Carlo algorithm. This constitutes the first example 
of a complete solution of the complex action problem 
coming from a topological term. First results from the 
dual simulation of the U(l) gauge-Higgs model with 
non-zero 9 were presented in [5]- 

The two-dimensional U(l) gauge-Higgs model is 
not only interesting as a testbed for the dual ap¬ 
proach, but also provides interesting physics. In its 


Higgs phase directly related to the phenomenologi¬ 
cal Ginzburg-Landau model of superconductivity [9 
it exhibits well-localized multi-vortex solutions with 
non-trivial integer topological charge ITOI - n^l . On the 
lattice the model and its topological effects have been 
studied in two |141 1 15] as well as in four dimensions 

mm- 

Here we are going to demonstrate that the topolog¬ 
ical features of the model in two dimensions can be 
further explored with the new techniques which now 
allow for simulations at arbitrary 8. 


II. U(l) GAUGE HIGGS MODEL WITH A 
TOPOLOGICAL TERM 


A. Continuum formulation 


The Euclidean continuum action for the 2- 
dimensional U(l) gauge Higgs model reads 


S[A, <j>\ = d 2 x (D ll (x)<j>(x))*D IJ ,(x)^(x) 


( 1 ) 


+ m 2 \(j){x)\ 2 + A|</>(:e)| 4 + ^-F tlv (x)F liV (x) 


Here x) = d^+iA^f x) denotes the usual covariant 
derivative and F^ v (x) = d tl A v (x)—d v A ll {x) is the field 
strength tensor for the U(l) gauge held A^x). The 
matter fields are described by a complex scalar held 

The topological charge for U(l) gauge helds in two 
dimensions is given by 

Q[A] = ^ Jd 2 xe^F^(x) = ^ Id 2 x F 12 (x) . (2) 


The topological term is added with a vacuum angle 8 




2 


to the action and the partition function reads 

Z = Jv[A]V[4>] e -S[AM-if>Q[A] _ (3) 

B. Lattice formulation 

On the lattice the partition function is given by 
Z = JV[U] V[cj>] e -s G [u]-s M [u^\-ieQ[u] _ 

For the gauge part we use the standard Wilson gauge 
action 


S g[U] = -f £[U*, P + U lv\ > (5) 

X 


where the sum runs over the sites a; of a 2-dimensional 
N t x N s lattice with periodic boundary conditions. 
The plaquette variable U X}P is built from the Un¬ 
valued link variables U Xtll , n = 1,2, and is given by 
U x , p = U x> i U x+ i 2 U* ~ 1 U * 2 ■ The matter part of 
the action reads 

Sm[U,<I>\ = [«I^| 2 +A|^x| 4 (6) 

X 

— ^^XPx^x^Px+p, + 4 > xU xfl (j) x+ p) , 

A* 


with a mass parameter k = 4 + m 2 . 

We stress at this point that in the literature (see, 
e.g., US) a different nomenclature can be found, 
where the couplings have a slightly different mean¬ 
ing (we add primes to distinguish them from our cou¬ 
plings): k' is the factor in front of the nearest neighbor 
terms, X' is the factor of a shifted quartic term of the 
form (10a:| 2 — l) 2 , and the quadratic term comes with 
a factor of 1. The connection between the two con¬ 
ventions is given by the transformations 


A = 


A' 
n' 2 


K = 


1 - 2X' 
k’ 



(7) 


and by dropping an irrelevant term. In the conven¬ 
tion [IS the limit A' —» oo freezes the radial mode to 
\<j)' x \ = 1 and one expects a Kosterlitz-Thouless tran¬ 
sition when varying k' . In our convention this cor¬ 
responds to following the line A = (1 — k'k)/2k' 2 for 
A —» oo (compare also the k-X phase diagram in Fig. [6] 
below). 

On the lattice the topological charge can be dis¬ 
cretized with, e.g., the “field-theoretical definition”, 

Q[U] = - K P ], ( 8 ) 


which reproduces ([2]) in the continuum limit. For the 
naive continuum limit this can be seen, by setting 
and noting that 

U xp = e l( ^- 1+A -+ i ^A I+ 2 ,i“ y W 2 ) = 1 + iF 12 (x) + ... 

(9) 

The measures T>[U\ and £>[</>] in the lattice path in¬ 
tegral are defined in the usual compact way, i.e., as 
product measures over all degrees of freedom on the 
lattice, 


Jm = n/' 


dA 


X,l± 


2t r 


S m = n / 0 


d(f)x 

2n 


( 10 ) 

Putting things together we can write the partition sum 
as 


Z = Jv[U]e v ^ + 71U Z* Z M [U\ , 

Z M [U\ = JV[<t>\ , (11) 


where Zm [U] is the partition sum of the matter fields 
in the gauge background. For a convenient notation 
of the terms that combine the gauge action and the 
topological charge we defined 


V = 


P 

2 


0 

47T ’ 



e 

47T 


( 12 ) 


It is obvious, that the conventional representation 


(11) of the lattice model is not suitable for a Monte 
Carlo simulation at 9 ^ 0, since then r] ^ rj and the 
Boltzmann factor is complex. In the next subsection 
we show that this problem is overcome by mapping 
the partition sum to dual variables. 


C. Dual representation 


By expanding the Boltzmann factors containing the 
nearest neighbor terms in Zm[U\ one can exactly map 
the partition sum of the matter fields into a dual form 
hm where the new degrees of freedom for the matter 
fields are loops dressed with the link variables U x , p . 
For the gauge fields, one proceeds in a similar way 
0 expanding the Boltzmann factor, rearranging 
terms and integrating out the link variables. A more 
detailed account of this exact transformation of the 
partition function Z into its dual form is provided in 
the appendix. 

In terms of the dual variables the partition function 
is given by 
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*= E 

{ij,p,p} 


n 


(|^rc,/i| H - ^,/i)* lx,n • 


n p( - 


n 


fj(\Px l+Px)/2+Pa: ^(|P*|-P*)/2+Px 


(|Px| + Px)! Pa:! 


^ f ^ [/x,/i /i,/i] ^ $(Px Px— 2 + ^x,l) £(Z> x _i p x ~\~ lx, 2 ) 


(13) 


In the dual representation the partition function is a 
sum over the set {l,l,p,p} of all configurations of the 
integer valued dual variables 


Px e 


lx, 11 5 Px F 0 , 


(14) 


x,/x and l x ^) or the 


which are assigned to the links (l 
plaquettes ( p x and p x ) of the lattice. For each config¬ 
uration there is a real and positive weight factor which 
consists of the terms in the first line of (13), where we 


have introduced the following abbreviations 


/•OO 

P{n x ) = dr r nx+1 e ~^ 2 ~ Xri , (15) 

Jo 

^x — ^ ( [ |^x,/a| T |^x— p,p\ ~b 2( l x ,/i 4~ l x —£ l,il ) ] ■ 

P 


rj < 0. In practice this is, however, an irrelevant re¬ 
gion of the parameters, since we are interested in the 
limit /3 —¥ 00 where one approaches the continuum 
limit (see below). Thus we can restrict ourselves to 
the parameter region /3 > 6/2tt where both r] and rj 
are positive. 

D. Observables and their dual representation 


For the analysis here, we focus on studying the be¬ 
havior of various bulk observables (for a calculation 
of propagators in the dual picture see, e.g., mm)- In 
particular we consider: 

Square of the absolute field and its susceptibility: 


' |2 > = 


-1 d 
N„ N t dn 


InZ 


X<t> 


1 d 2 

N s N t d^ 


Plaquette and plaquette susceptibility: 


InZ . 
(16) 


A subset of the dual variables, i.e., the l X:/J and the 
p x , are subject to constraints, which are collected in 
the second line of (13). Here the 5 denote Kronecker 
deltas, i.e., S(n) = 0 ^, 0 - The constraints for the l x )(U 


are a discretized version of =0 and thus imply the 
conservation of /-flux at every site x. The remaining 
constraints are associated with the fluxes along the 
links of the lattice: Here the flux at a link introduced 
by a non-trivial p x has to be compensated by an op¬ 
positely oriented flux from a neighboring plaquette or 
by /-flux from link variables l x The combinations of 
all constraints gives rise to admissible configurations 
that contribute to the partition sum, which consist 
of closed loops of /-flux which are filled with occu¬ 
pied plaquettes such that at each link the total flux 
is zero. These admissible configurations are a natu¬ 
ral reduction of the admissible configurations that are 
discussed in more detail for the 4-dimensional case in 


m\M- 


We close this subsection with remarking that the 
weight factors in the dual representation (13) are al¬ 
ways real, but become negative when either r) < 0 or 


(Re U p } = 


1 d 
N s N t dp 


d 2 


lnZ ’ Xp 55 N^N t W XnZ - 

(17) 

Topological charge density and topological charge 
susceptibility: 

— 1 Q —\d 2 

(</> = > Xt = TTTFwxyln^. (18) 


N s N t 86 


Ns N t d9 2 


The dual expressions of these observables can be ob¬ 
tained by evaluating the derivatives of In Z using the 
dual representation of the partition function Z . We 
give two examples for dual expressions of observables, 
the somewhat simpler field expectation value, 


(H 2 > 


d _ —1 /y- dP(n x ) 1 \ 

N s N t dn n N s N t \Z^ 8 k P{n x ) J 


1 

N a N t 



P{n x + 2)\ 
P{n x ) j ’ 


(19) 


and the slightly more involved expression for the topo¬ 
logical charge susceptibility, 


Xt = 


-1 

NsN t 


(47 nf) 2 






1 / 


8n 2 r]r] \ 




M -s + s 

2 2 


( 20 ) 


(47T77) 2 


[Sl-l+sl 

[!|!-f + 3-il 

W-T 

[If-f+sl 

1 

[H+f+sl 

yi 

2 2 

2 2 

/ \ 47t?7 

2 2 

47t?7 

2 2 

/ J 
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where we use the abbreviations |<S| = \p x \, S = 

'^2 x Px and S = J2xPx f° r various sums of the plaque- 
tte occupation numbers p x and p x . The expectation 
values on the right-hand sides of (19) and (20) are un¬ 
derstood as expectation values in the dual represen¬ 
tation. In a similar way all bulk observables defined 
above can be obtained as weighted moments of the 
dual variables. 


III. TESTS IN PURE GAUGE THEORY AND 
DUAL MONTE CARLO UPDATES 

A. Pure U(l) gauge theory: Semi-analytical 
results and continuum limit 


If we neglect the matter fields (’’quenched case”), 
the model can be solved (semi-) analytically, i.e., we 
obtain for the partition function a simple and fast con¬ 
verging sum, which can be evaluated efficiently to ar¬ 
bitrary precision. 

For the pure gauge system with topological term 
the partition function reduces to 


* = £ 
M 


nz 


(Vvv) 


=\P* + |Px| 


x p^O +Px ) ! Px[ 



X 


II S (.Px-Px-2) S (P X -i - Px) , 


( 21 ) 


where the sums over the p x in the first parentheses 
are well known and yield the modified Bessel func¬ 
tions I n (x). We thus obtain for the partition sum the 
expression 


* = £ 

fp> 


n j ip*i( 2 ^) 



X 


II S (Px-Px-2) S (p x -i -Px) 


it numerically with Mathematica and the results we 
show below were obtained in this way. 

Eq. (231 nicely illustrates how the vacuum angle 9 


influences the physics in our sy stem . For 9 < 0 we 
have r] > rj and thus the term ( y/p/fj) q enhances con¬ 
figurations with q > 0. Configurations with q x = q\/x 
correspond to configurations of constant electric flux 
and via the term (y/p/fj) 9 the vacuum angle allows 
to introduce such flux in the system. 

We can also use (23) for a first assessment of the 


continuum limit. In particular it is interesting to 
study how well the field theoretical lattice definition 
of the topological charge ^ can reproduce the ex¬ 
pected 27r-periodicity of observables as a function of 
the vacuum angle 9. This is not a priori clear, since 
the definition ([8]) does not guarantee an integer val¬ 
ued topological charge - this is expected only in the 
continuum limit. The continuum limit is approached 
via 


/3 —► oo iV s , N t —► oo with 


P 


N s N t 


= const. (24) 


Dimensional analysis yields \B\ = L 2 which implies 
that the continuum limit \2A\ corresponds to keep¬ 
ing a fixed physical volume. One expects that in the 
fixed-volume continuum limit observables become 2-7T- 
periodic in 9. 

We studied the 0-dependence of various observ¬ 
ables and as an example in Fig. |T] we show the pla- 
quette expectation value as a function of 9 on a se¬ 
quence of lattices that approach the continuum limit 
(/? = 0.1,2.5,10.0 and 40.0 at fixed /3/N s N t = 0.001). 
The results were obtained by evaluating (ReU p ) = 
In Z/N t N s for the pure gauge partition sum as 


given in (23). The tests documented in Fig. [I] show 
that indeed 27r-periodicity is recovered and the field 
theoretical definition ([8]) is well suited to study the 
continuum limit of the model with the topological 
term. 


B. Dual Monte Carlo simulation 


In the quenched case we have no matter flux for sat¬ 
urating the constraints at the links. Thus the Kro- 
necker deltas in (22) force the plaquette occupation 
numbers p x £ Z to have the same value q at each 
lattice site such that all fluxes along the links cancel. 
Hence every configuration that obeys all constraints 
can be labeled by a single integer q and p x = q \/x. 
This simplifies the partition sum to 


Z = 


+oo 


q =—oo 


7 kl( 2 Vw) 


N,N t 


(23) 


The modified Bessel functions I n decay faster than 
exponentially with the index n and thus the series (231 


converges rapidly. It is straightforward to evaluate 


Before we come to the presentation of the results 
for the full model, let us discuss the Monte Carlo up¬ 
date strategy used for the dual representation. In the 
dual representation the dynamical degrees of free¬ 
dom are the integer valued plaquette and matter flux 
variables, p x ,p x and l x ,p,,lx,p■ The dual variables p x 
and l x ^ are not subject to any constraints and they 
can be updated independently in the usual way with 
a local Monte Carlo update scheme. More demand¬ 
ing are the variables p x and l XjP- which have to obey 
the constraints of conserved /-flux at each site and a 
vanishing of the combined p- and /-flux at each link of 
the lattice. 

Although it is possible to find a generalization of 
the worm strategy m to abelian gauge-Higgs models 
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[3 = 0.1, N s =N t =10 [3 = 2.5, N s =N t =50 




0 


0 


[3=10, N s =N t =100 



-4n -3n -2n -In On In 2n 3n 4n 
0 


(3=40, N s =N t =200 
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FIG. 1. Plaquette expectation value (Re U p ) of the pure gauge theory versus the vacuum angle 8. We show the approach 
to the continuum limit using /3 = 0.1, 2.5,10.0 and 40.0 at fixed /3/N s Nt = 0.001. The plots nicely demonstrate how the 
observable becomes 27r-periodic in 8 when approaching the continuum limit. 


in arbitrary dimensions 012, for the 2-dimensional 
model studied here we use a simpler local update for 
the dual variables. The update strategy consists of 
two types of updates: 

1. A local plaquette/link update: For a lattice site 
x we randomly choose A x = ±1 and propose to 
change 


Px -t p' x = Px + A x , 

lx,\ l x ,l ~ A X , ^x+1,2 l x+i,2 _ > 

^ x + 2,1 ^ C + 2,1 + 5 ^,2 C ,2 + ^ x ■ 

The change is accepted with a Metropolis step. 

2. A global winding gauge update: A A x = ±1 is 
chosen randomly and we propose to change 

Px -> p' x = Px + A Vx . (27) 


(25) 

(26) 


Again the change is accepted with a Metropolis 
step. 


It is easy to see that these updates leave the con¬ 
straints intact and that this procedure is ergodic. To 
be precise, the plaquette update alone is already er¬ 
godic, but mixing sweeps of the local plaquette/link 
update with global winding gauge updates consider¬ 
ably reduces the auto-correlation of observables re¬ 
lated to the topological charge. For the results we 
show here, we typically use 5 x 10 4 such combined 
sweeps for equilibration followed by 10 6 measurements 
separated by 5 combined sweeps for decorrelation. 
The error bars we show are statistical errors deter¬ 
mined with a jackknife analysis. 

The performance of the Monte Carlo updates (251, 
(f26| ) was studied for the 4-d case in great detail in 
( 2 , and qualitatively the performance behavior is the 
same in the 2-d case, which is, however, numerically 
considerably less demanding. As remarked, the up¬ 
date (27) is not necessary for ergodicity, but helps 
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P = 10, N s =N t =10 



0 

k=io i— •— i pure gauge 


FIG. 2. Plaquette expectation value (Re U p ) of the full 
model versus the vacuum angle 9. We show data for a large 
value of k = 10.0 at A = 1 N a = Nt = 10. We compare 
the numerical data from the Monte Carlo simulation to the 
results from pure gauge theory discussed in the previous 
section. 


for decorrelating topological quantities since it corre¬ 
sponds to changing the topological sector. Thus the 
acceptance of this update depends exponentially on 
the volume, which is, however, not surprising, since 
the action for an additional charge is extensive for the 
2-d U(l) case and the exponential volume-dependence 
is physical (this is different for SU(3) in 4-d). 


IV. NUMERICAL RESULTS FOR THE FULL 
MODEL 

A. Continuum limit and periodicity in the 9 
angle 


continuum limit and to compare the results for three 
different values of the mass parameter, n = 4.0, k = 
5.0 and again the heavy mass value k = 10.0. As in the 
case of pure gauge theory we consider the continuum 
limit for /3 —► oo with /3/N s N t fixed, and use /3 = 1.6, 
3.6, 6.4 and 10.0 at fixed /3/N s N t = 0.1. In Fig. [3]we 
show the results for the plaquette expectation value as 
a function of 9. As before in the quenched case also 
here we observe that, as we approach the continuum 
limit, the 27r-periodicity of the observable emerges as 
expected. It is remarkable that for this pure gauge 
observable the dependence on the mass parameter n 
is only very weak and the behavior of the plaquette 
in the full theory is nearly the same as for the pure 
gauge case as long as we are in the symmetric phase. 

Finally, in Fig. [4] we study the expectation value 
of the topological charge density ( q ) as a function of 
0, again using the values k = 4.0, 5.0 and 10.0 and 
approaching the continuum limit with the sequence 
/ 3 = 1.6, 3.6, 6.4 and 10.0 at fixed j3/N s N t = 0.1. As 
before we observe the emergence of 2-7r-periodicity in 
the continuum limit. 

Here we remark that the topological charge density 
(q) is a variable that is odd in 9, while the plaque¬ 
tte (Re U p ) is an even function. This can be seen 
from the fact that when changing the link variables 
in the path integral according to U XtP -A U* p we 
have Q\U] -A —Q\XJ], while R eU p remains unchanged 
(compare Eqs. 0 and ®). This leads to a lin¬ 
early rising behavior of (q) in the vicinity of 9 = 0. 
This (anti-) symmetry is clearly visible in all plots of 
Fig-0 while the expected 27r-periodicity in 9 is recov¬ 
ered only in the continuum limit. 

As already for the plaquette expectation value, also 
for the topological charge density (q) we find that the 
results are essentially independent of the mass param¬ 
eter k. This indicates that the physics of the model 
related to the vacuum angle 9 is dominated by the 
behavior of the gauge sector and the dynamics of the 
matter field plays only a sub-leading r ole as long as 
we are in the symmetric phase (c.f. Sec. IV B). As we 
will see below, the behavior changes considerably in 
the Higgs phase. 


Now we study the full model with matter fields. As 
a first consistency check we discuss a simulation at 
large k, i.e., large mass, where the matter fields be¬ 
come static and the results are expected to approach 
the quenched results from the last section (the quar- 
tic coupling A is always set to A = 1 in this sub¬ 
section). Figure [2] shows the plaquette expectation 
value as a function of 9 for k = 10.0, /3 = 10.0 and 
N s = N t = 10. It is obvious that, as expected, the 
numerical data from the dual simulation are indeed 
in very good agreement with the curve obtained from 
explicitly summing the quenched partition function 

Now we switch to lighter matter fields to study the 


B. Phase diagram in the k-A plane 

We now discuss the determination of the phase di¬ 
agram in the plane of the mass parameter n and the 
quartic coupling A. In more than two dimensions 
one expects a true phase transition line separating 
the symmetric from the Higgs phase. In two dimen¬ 
sions one only finds smooth behavior of the variables 
since spontaneous breaking of a continuous symmetry 
is not possible in two dimensions due to the Mermin- 
Wagner-Coleman theorem. Second order derivatives 
of the free energy show an extremum but no scaling 
with the volume. 
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FIG. 3. Plaquette expectation value (Re U p ) of the full model versus 0 for three different values of the mass parameter 
k at A = 1. We show the approach to the continuum limit using (5 = 1.6, 3.6, 6.4 and 10.0 at fixed /3/N a Nt = 0.1. The 
discontinuity in the top left plot near 0 = ±3-7r is due to the violation of (5 > 9/2-k. 


An example is given in Fig. [5] where we show the 
field expectation value ( \<j >\ 2 ) (lhs. plot) and the cor¬ 
responding susceptibility \4> (rhs.) versus the mass 
parameter n at A = 0.1 and /3 = 10.0 for different vol¬ 
umes. It is obvious that the data for different volumes 
fall on top of each other and volume scaling is absent. 

Nevertheless, one expects physical properties to 
change when crossing the line of the smooth transition 
and in order to determine the location of this line in 
the k-X plane, i.e., the phase diagram, we identify the 
position of the peak in higher derivatives for various 
values of the parameters. The results are shown in 
Fig. [G} For the determination of the transition line 
in Fig. [6] we used the maxima of x<t> as a function 
of k (plusses) and the maxima of d/dn(ReU p } as a 
function of both, k (crosses) and of A (squares). We 


stress, that at a crossover transition different higher 
order derivatives do not necessarily have to peak at 
the same position. However, as is obvious from Fig. [6] 
we here find agreement among all observables within 
error bars. We find a line of crossover points extend¬ 
ing from k = 4.0 and A = 0 (which corresponds to the 
free massless field) to n = 2.0 and A = 1 (the largest 
value of A considered here). The data for Fig. [6] are 
for N s = N t = 10 and /3 = 10.0, i.e., relatively close 
to the continuum limit. 

An interesting question is how strong the position 
of the crossover line depends on the vacuum angle 9. 
In our study we find that this dependence is only very 
weak. This is illustrated in Fig. [TJ where the peak 
of x\ 4 > I 2 is shown in high resolution for four different 
values of 9 G [0,27r] at A = 0.5. Starting at 9 = 0 the 
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FIG. 4. The topological charge density (q) versus the vacuum angle 0 for three different values of the mass parameter k 
at A = 1. We show the approach to the continuum limit using (3 = 1.6, 3.6, 6.4 and 10.0 at fixed /3/N s Nt = 0.1. Again 
the discontinuity in the top left plot near 9 = ±37r is due to the violation of /3 > 8/2n. 


peak gets shifted towards slightly larger values of the 
mass parameter (here from n ss 2.84 to k ~ 2.88) until 
9 = 7r, after which, due to periodicity in 9 , the peak 
moves back to its position at 9 = 0 when we reach 
9 = 2n. The situation is similar for other values of the 
coupling A € [0,1]. We conclude that the variation of 
9 amounts to only a very small shift of the crossover 
line. 


C. Characterization of the two phases 

As a next step in our explorative study of using 
dual variables in the U(l) gauge-Higgs model with 
topological term we attempt a characterization of the 
two phases - in particular with respect to their de¬ 
pendence on the vacuum angle 9. We remark again 
that the transition is smooth (it is of the Kosterlitz- 


Thouless type) and no simple order parameter exists. 
Nevertheless, for completeness, we start with summa¬ 
rizing the characteristic behavior of all observables in 
the two phases for 9 = 0. The corresponding results 
are presented in Fig. [8j where the six observables 
discussed above are shown as a function of the mass 
parameter n for A = 0.5. In other words, we analyze 
horizontal slices through the phase diagram Fig. [6] and 
we compare the results for four different volumes. In 
all observables a changing behavior can be seen near 
n = 2.8, the value on the critical line for A = 0.5. 
The topological charge is (within error bars) identi¬ 
cally zero, as expected at 9 = 0, but obviously the 
fluctuations are different in the two phases, as can be 
seen by the change of the size of the error bars and the 
increase of the fluctuations. The observables do not 
depend on the lattice volume, except for \top which 
shows a small finite volume effect in the symmetric 
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FIG. 5. Expectation value of the square of the field modulus (|<^>| 2 ) and the corresponding susceptibility versus the 
mass parameter /t at A = 0.1 and ft = 10.0 for different volumes. One can clearly see the changing behavior of the first 
derivative starting at approximately k = 3.7, leading to a peak in the corresponding susceptibility. The observable is 
essentially independent of the volume indicating a smooth transition. 



K 


FIG. 6. Phase diagram in the plane of couplings A and 
k. We find two phases separated by a crossover line. The 
phase boundary was determined from the maxima of x<!> 
as a function of k (plusses connected with a dotted line) 
and the maxima of 9/<9 k(R eU p ) as a function of both, k 
(crosses) and of A (squares). 


phase for the smallest lattice size (see also the dis¬ 
cussion below). Let us notice at this point that the 
behavior of the topological susceptibility versus k for 
fixed A and /3 at 6 = 0 has been checked by simulat¬ 
ing the model with the standard Metropolis algorithm 
without dualization of the variables. We found exact 
agreement. 

At this point we remark that the topological sus¬ 
ceptibility at 9 ^ 0 can also be negative, despite the 



K 


FIG. 7. Variation of the peak in \<t> with 6 at A = 0.5, 
P = 10.0 and N s = N t = 10. The vertical lines delimit the 
range of K-values where the cross-over is found (as deter¬ 
mined from the maxima of \<t >) when varying 6 between 
6 = 0 (full red line) to 6 = 7r (dashed blue line) and back 
to 6 = 2n. 


fact that it is a second derivative of In Z. This can be 
seen as follows: Instead of integrating in the path inte¬ 
gral over the variables we can also integrate 

over the complex conjugate variables U* jt , (f>*. Us¬ 
ing the symmetry properties (compare equations @ 






















10 






FIG. 8. Observables versus the mass parameter k at A = 0.5, /3 = 10.0 and 9 = 0 for different volumes. One can clearly 
see the changing behavior of all observables at approximately k = 2.8 corresponding to the crossover position at A = 0.5 
as mapped out in the phase diagram in the previous section. 


- §) JV[U*} = JV[U], fV[<f> *] = JV[4>], Sg[U *] = 
sm, Sm[U*,P} = S M [U, (j)\ and Q[U*] = -Q[U ] 
we find for the partition sum 

Z = Jv[U\D[<P] e -SG[u]-s M [u,<t>\ e -i8Q[u] ^8) 

= J V[U]V[cj)} e ~ So[V)~ Sm[u,4>\ C0S (6»Q[t/]) . 

Evaluating the topological susceptibility as the second 


derivative of In Z with respect to 9 (compare Eq. ([18])) 
we obtain 

N s N t Xt = -g 0 . 2 hi Z (29) 

= ^ JV[U]V[<t>] e - Sg{u]-Sm[u,<i>] cos (9Q[U]) Q[U} 2 

+ (^ Jv[U]V[^}e- s ^ u] - SM ^sm(9Q[U})Q[U}^ . 
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The right hand side has two contributions, with the 
second term being the square of the first derivative 
and thus positive. The first term, however, comes 
with cos(0Q[[/]), i.e., the even part of exp (—id Q[U]). 
This first term can be negative, and thus also Xt needs 
not be positive. The same of course also holds for the 
dual expression of Xt given in (201. 

In our analysis of the model with the dual variables 
approach we find that the two phases can be charac¬ 
terized by their response to the variation of the vac¬ 
uum angle. As a first illustration of this fact in Fig. 
[9] we show 3-D plots of Q and Xt as a function of k 
and 9 at A = 0.5, i.e., when considered as a function 
of k we again inspect a horizontal slice through the 
phase diagram Fig. [6] and we expect to see changing 
behavior at k ~ 2.8. This is indeed what we observe: 
At large values of the mass parameter n, i.e., in the 
symmetric phase, we see oscillatory behavior with 9 in 
both observables, while in the crossover region the ob¬ 
servables are independent of 6 within error bars. The 
transition between the two types of behavior takes 
place as expected near k ~ 2.8. 

In Fig. 10 we now look at the volume dependence 
of Xt at 6 = 0 and 9 = tt and compare the behavior for 
different values of the mass parameter k. The topo¬ 
logical charge density (q) has a negative slope at 6 = 0 
and therefore Xt is negative on the lhs. plot. We ob¬ 
serve a strong dependence on the mass parameter k. 
The susceptibility essentially vanishes in the broken 
phase and then starts to deviate from 0 for k > 2.5 . 
For all values of k a saturation is reached on lattice 


volumes between N s = N t = 10 and N s = N t = 12. 
For 9 = 7r (rhs. plot) the behavior is different: Here Xt 
also vanishes in the broken phase and then is positive 
for k > 2.5. Most remarkably, Xt does not seem to 
reach saturation as a function of the volume, a fact 
that hints at a possible phase transition. 


D. The transition at 9 = it 


In the previous section we found evidence that in 
the symmetric phase there might be a transition at 
9 = tt. To identify the transition we analyzed the In¬ 
dependence of (q) and Xt for a point in the symmetric 
phase, in particular at A = 0.5, k = 4.0 with /3 = 10.0. 
The results as a function of 9 for different volumes are 


shown in Fig. 11 It is obvious that xt (rhs. plot) has 


a maximum near 9 = tt, and the height of the peak at 
the maximum increases with the volume. A detailed 
analysis shows, that the height of the maximum of Xt 
scales almost perfectly with the volume, which indi¬ 
cates a first order transition. This is reflected in the 
behavior of (q), which in the large volume limit de¬ 
velops a discontinuity near 9 = tt. The analysis was 
repeated at other points in the symmetric phase with 
the same result and we conclude, that in the symmet¬ 
ric phase the system has a first order phase transition 


as a function of 9. 

The transition at 9 = tt can be related to charge 
conjugation, i.e., the discrete symmetry transforma¬ 
tion 17 X>(U —> U a * Vx,fi and <j) x -X (j>* Wx. While the 
action is invariant under this symmetry, the topologi¬ 
cal charge changes sign: Q[U] —> —Q[U}. However, 
at 9 = n the Boltzmann factor with the topologi¬ 
cal charge is given by = (— l)^!^] which is 

again symmetric under Q[U] —> —Q[U\ and thus un¬ 
der charge conjugation. This discrete symmetry can 
be broken in the symmetric phase of the model. 

A careful inspection of Fig.[ll]reveals that the tran¬ 
sition is not located at exactly 9 = tt, but appears at 
a slightly larger value of 9. This can be attributed 
to an effect of finite lattice spacing. To better under¬ 
stand this behavior we looked at the observables (q), 
Xt, also in the pure gauge case where we can use the 
semi-analytical results from Section III A The corre¬ 
sponding results for /3 = 2.5 and for /3 = 10.0 on four 
different volumes are shown in Fig. |12[ For the value 
/3 = 10.0 (plots in the bottom of Fig. |12|) we observe a 
behavior very similar to the one of Fig. 11 However, 
at (3 = 2.5, where we are further away from the con¬ 
tinuum limit, the transition is seen at an even larger 
value of 9 and indicates that the position 9 cr n = tt for 
the transition is reached only in the continuum limit, 
i.e., for /3 —► 00 . 

Only in the continuum limit Q[U] becomes re¬ 
stricted to integers, while away from the continuum 
limit Q\U] has a distribution that is localized not ex¬ 
actly around integers, but has its maxima at values 
shifted to slightly smaller numbers than the integers 
of the continuum limit. For example in the charge 1 
sector one could find the maximum at Q ma x{P) = 0-8 
instead of the continuum limit value Q m ax{o o) = 1. 
The shift of the transition towards values of 9 larger 
than tt then is explained by the condition that 9 cr i t x 
Qmax(P) = TT, such that the symmetry U x ^ -A U*^ 
emerges. This gives rise to 9 cr u = tt /Qmax{/3), which 
explains 9 cr n > tt for B < oc. 

We remark, that we also looked at other observables 
in the pure gauge case, in particular (Re U XJ) ) and Xp- 
Also there we find a first order behavior at 9 = tt. 


V. SUMMARY AND DISCUSSION 

In this exploratory study we explore strategies for 
using dual variables in a simulation of a lattice field 
theory with a topological term. We show for the 
case of the U(l) gauge-Higgs system in two dimen¬ 
sions (scalar Schwinger model) that a dual formula¬ 
tion can be found where the complex action problem 
of the conventional representation is overcome. The 
dual variables are loops for the matter fields and sur¬ 
faces for the gauge fields and the partition sum has 
only real and positive contributions. We show that 
in terms of the dual variables a Monte Carlo Simula- 
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FIG. 9. Topological charge density (Ills, plot) and topological susceptibility (rhs.) versus k and 6 at A = 0.5, fi = 10 and 
N s = Nt = 10. 
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FIG. 10. Topological susceptibility versus the volume V = N s x Nt for different k at A = 0.5 and fi = 10 at 8 = 0 (lhs) 
and 8 = it (rhs). 


tion is possible at finite values of the vacuum angle 9. 
This constitutes the first example of a theory with a 
vacuum angle where a simulation could be performed 
with a complete solution of the corresponding complex 
action problem. 

Using the dual approach we show that for the 
plaquette expectation value and for the topological 
charge (Q) the expected 27r-periodicity emerges in 
a properly implemented continuum limit. This was 
shown for the semi-analytically tractable case of pure 
gauge theory, as well as for the full theory with mat¬ 
ter fields. For the latter case we also found that the 
behavior of the observables as a function of 9 is essen¬ 
tially independent of the mass parameter k in the sym¬ 


metric phase. In the broken phase we find a quantita¬ 
tive dependence on the mass parameter but no qual¬ 
itative one, i.e., observables acquire a constant shift 
on top of which they show exactly the same periodic 
behavior as in the unbroken phase. A clear distinction 
between the Higgs- and the symmetric phase appears 
for 9 = 7r, where we identify a first order behavior of 
observables in the symmetric phase, which is absent 
in the Higgs phase. 

In |17] the same model was studied, however in four 
dimensions with an external source. There it was 
found that the Higgs phase splits into two regions, dis¬ 
criminated by the kind of magnetic flux penetration. 
Therefore we covered a large range of mass param- 
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FIG. 11. Topological charge and xt versus 6 at k = 4.0, A = 0.5 and /3 = 10.0 for different lattice volumes. 






FIG. 12. ( q) (lhs.) and (rhs.) in pure gauge theory as a function of 9. We show results at /3 = 2.5 (top) and (5 = 10.0 
(bottom) and compare four different volumes. 
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eter values k £ [—50,4] for A = 0.5 to search for a sum can be directly taken over from [TH1] 
hint of a phase change in the broken phase. However, 
the only transition behavior we saw, was the crossover 
from the symmetric to the Higgs phase as discussed 
in Section |TVB[ The topological charge as well as the 
topological susceptibility stay constant (within error 
bars) throughout the broken phase. 


Z M [U] = Y, 
{i,T} 


n 


(]4,ft| + 4,y)! 4 iAt !) 




n^) 

X 

, (30) 


APPENDIX 


Derivation of the dual representation 


with l x u £ Z, l x fi £ INq and P{n x ) and n x as defined 
in Eq. (151. In its dual form Zm[U] is a sum over loops 
of conserved 1-flux and the loops are dressed with the 
link variables via the terms (L4 iA1 ) . 

The full partition sum Z is obtained by integrating 
Zm[U] over the gauge fields U Xtft with the Boltzmann 
factor e~ SG ^ u \ The dressed loops in (301 provide the 
additional gauge field dependent factor Tl,,, 
such that we need to compute the integral 


Z G [l\ = Jv[U} 

= [ V[U] 


X,[A 

n w x j‘ 


o~S G [U]—i 0 Q[U) 


(31) 


n 


,v u x , p e vu: 


In this appendix we provide a brief summary of the 
steps leading to the dual representation (13) of the 
partition sum. Initially the derivation as given for 
abelian scalar and gauge fields in four dimensions in 
but the adaption of the arguments to two di¬ 
mensions and the topological term is straightforward. with rj and fj defined in (12). The dualization of the 

gauge fields proceeds in a way equivalent to the mat¬ 
in (11) we have defined the partition sum Zm[U] ter fields: The Boltzmann next step is to expand the 

of the matter fields in a background configuration of exponentials at each lattice site and for the two terms 

lattice gauge fields U X:fl . For Zm[U] the dual partition in the exponent separately 


Z G [l\= V[U] 


OO OO 

n<«W"n E E h 


x,n 

OO 


x n x = 0 


n x — 0 


n e 

_ x n x ,n x = 0 . 


n 


ryjTl'X jjTl x 

n x \ n x \ 


V[U) 


^ rJit-n*— n x _2+n x _$+l x ,i (jj ^ n a: —n a: —n x _f+n a ._i+l x ,2 


= E 

{ n,n } 

= E 

{ n,n } 


n 

X 

n 


n x \ h x \ 


j-j^x y~j^ x 

n x \ n x \ 


n 


^dA x>1 




T dA 


2i r 


x ’ 2 c iA x .2( n x -n x -n —+n^_, +^.2) 


2tt 


S ( n x - n x - n x _2 + n x _% + 4,i) S(n x - n x - n x _i + n x _i + 4, 2 ) 


(32) 


r 


Upon defining new summation variables p x £ Z and final form 
Px € INq 


Px — Tl x Tl x 3 \Px\ + 2 p x — Tl x Tl x , ( 33 ) 


z g[1\ = E 

{p,p} 


n 


? y(|P^I+Px)/2+p a: 7jj(\pi:\-pA/ 2 +Px 


(\Px\ +Px) ] - Px ! 

}js(px ~Px- 2 +4,l) s(p x _ 1 -p x + 4,2 


(34) 


we get the partition function for the gauge fields in its Combining this result with the matter part from (30) 















































completes the derivation of the dual representation 
and yields the expression of Eq. (13). 
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